Fischer Decomposition for Difference Dirac 

Operators 

N. Faustino* 
Departamento de Matematica, 
Universidade de Aveiro, 
P-3810-193 Aveiro, 
Portugal 

U. Kahlert 
Departamento de Matematica, 
Universidade de Aveiro, 
P-3810-193 Aveiro, 
Portugal 

February 2, 2008 



Abstract 

We establish the basis of a discrete function theory start- 
ing with a Fischer decomposition for difference Dirac opera- 
tors. Discrete versions of homogeneous polynomials, Euler and 
Gamma operators are obtained. As a consequence we obtain a 
Fischer decomposition for the discrete Laplacian. 
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1 Introduction 



Clifford analysis is a powerful tool to solve some kinds of problems related with 
vector field analysis. 

A comprehensive description of Clifford function theory was introduced by F. 
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Brackx, R. Delanghe and F. Sommen in ^ and later by R. Delanghe, F. Sommen 
and V. Soucek in p]. 

In [HEI, K. Giirlebeck and W. Sproi^ig proposed some strategies to solve 
boundary value problems based on the study of existence, uniqueness, repre- 
sentation, and regularity of solutions with the help of an operator calculus. In 
the same books, the authors introduce also the basic ideas to develop a discrete 
counterpart of the continuous treatment of boundary value problems with the 
introduction of a discrete operator calculus in order to find a well-adapted nu- 
merical approach. An explicit discrete version of the Borel-Pompeiu formula 
was presented for dimension n = 3. 

This was further developed in where K. Giirlebeck and A. Hommel 

developed finite difference potential methods in lattice domains based on the 
concept of discrete fundamental solutions for the difference Dirac operator which 
generalizes the work developed by Ryabenkij in JHI- A numerical application of 
this theory was presented recently by N. Faustino, K. Giirlebeck, A. Hommel, 
and U. Kahler in P] for the incompressible stationary Navier-Stokes equations. 
In this paper, the authors proposed a scheme which solves efficiently problems 
in unbounded domains and show the convergence of the numerical scheme for 
functions with Holder regularity which is a better gain compared with the con- 
vergence results for classical difference schemes. 

Moreover, while all these papers claim to be based on discrete function theo- 
retical approaches, from the concepts of the theory of monogenic functions only 
the Borel-Pompeiu formula and with it Cauchy's integral formula were obtained. 
There is no "real" development of a discrete monogenic function theory up to 
now. 

This paper is supposed to be a step in this direction. To this end discrete 
versions of a Fischer decomposition, Euler and Gamma operators are obtained. 
For the sake of simplicity we consider in the first part only Dirac operators which 
contain only forward or backward finite differences. Of course, these Dirac 
operators do not factorize the classic discrete Laplacian. Therefore, we will 
consider in the last chapter a different definition of a difference Dirac operator 
in the quaternionic case (c.f. [3) which do factorizes the discrete Laplacian. 

Let us emphasize in the end a major obstacle in the discrete case. While in 
the continuous case the are only one partial derivative for each coordinate Xj 
we have two finite differences in the discrete case. Therefore, we will have not 
only one Euler or Gamma operator as in the continuous case, but several. Each 
one will turn out to be connected to one particular Dirac operator. 

2 Preliminaries 

Let ei, . . . , e„ be an orthonormal basis of R". The Clifford algebra CTo,n is the 
free algebra over M" generated modulo the relation 

= -|a;peo, 
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where eo is the identity of C£o,n- For the algebra C£o,n we have the anti- 
commutation relationship 

where Sij is the Kronecker symbol. In the following we will identify the Euclidean 
space K" with /\^ CTo.n, the space of all vectors of C£o,n- This means that each 
element x of M" may be represented by 

n 

X ^ ^ X^Si . 
i=l 

From an analysis viewpoint one extremely crucial property of the algebra C£o,n 
is that each non-zero vector a; G M" has a multiplicative inverse given by . Up 
to a sign this inverse corresponds to the Kelvin inverse of a vector in Euclidean 
space. Moreover, given a general Clifford number a = e^a^i, ^ C {1, . . . , n} 
we denote by Sea = a0 the scalar part and by Veca = eiai -I- ... -I- e„a„ the 
vector part. 

For all what follows let fl C M" be a bounded domain with a sufficiently 
smooth boundary F — dfl. Then any function f : fl i-^ C^o.n has a represen- 
tation / = J^A^^fA with M-valued components Ja- We now introduce the 
Dirac operator D — X)"^! T^^^s operator is a hypercomplex analogue to 

the complex Cauchy-Riemann operator. In particular we have that — —A, 
where A is the Laplacian over E". A function f : Q >~> C£o,n is said to be 
left-monogenic if it satisfies the equation {Df){x) — for each x G ft. A similar 
definition can be given for right-monogenic functions. Basic properties of the 
Dirac operator and left-monogenic functions can be found in Q], [5], [E], and j^. 

Now, we need some more facts for our discrete setting. To discretize point- 
wise the partial derivatives ^ in the equidistant lattice with mesh width h > 0, 
RJJ — {m,h = (mih, . . . ,mnh) : m £ Z"}, we introduce forward/backward 



differences 9^ 



n±i I ^,^ u{mh) - u{mh ± hei) 

dt'u{mh) = (1) 

These forward/backward differences d^^ satisfy the following product rules 

(a±7<?)(m/i) = f{mh){d^'g){mh) + (d^'f){mh)g{mh±he,), (2) 
{d^'fg){mh) = f{mh±he,){d^^g){mh) + {d^'f){mK)g{mh). (3) 

The forward/backward discretizations of the Dirac operator are given by 

n 

Dt^Y.^.dr. (4) 

1=1 

In this paper we will also use the following multi- index abbreviations: 
(m/i)(") := (mi/i)"i (ma/i)"^ . . . (to„/i)""; 
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a! := a'i!a2! . . . a„!; 
\a\ := ai + a2 + ■ . .an 

for a = (ai,Q!2, • • • = ejO!,. 

3 Fischer decomposition 

The basic idea of a Fischer decomposition is to decompose any homogeneous 
polynomial into monogenic homogeneous polynomials of lower degrees. In the 
classic case such a decomposition is based on the fact that the powers are 
homogeneous and that = sx^~^ . A first idea would be to consider instead 
of simply the powers {mhy , but while these powers are still homogeneous 
the last condition is not true in the discrete case, unfortunately. Therefore, 
we will start by introducing discrete homogeneous powers which will play the 
equivalent role of in the discrete case. 

3.1 Multi-index factorial powers 

Starting from the one-dimensional factorial powers 

s-l 

{rriih)'^'^ := 1, {rmh)^:^^ := {mih Tkh),seN (5) 

fc=0 

we introduce the multi-index factorial powers of degree |a:| by 

n 

{mhp=l[{mihp\ 

The one-dimensional factorial powers [rriih)^^ have the following properties 
PI. (mi/i)^p'+^^ = {mihTsh){mih)%^- 
P2. d^^{m^h)%^ = s{m,h)^^-^^5^j; 
P3. d^'irmh)';^^ =s{mihTh)''4~^^Sij; 

P4. (mi/i)^^ -^xt = {rriihy for h ^ 0, 

where 6ij denotes the standard Kronecker symbol. 

As a direct consequence of these properties, we obtain the following lemmas: 
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Lemma 3.1 The multi-index factorial powers of degree \a\, {mh)^^\ satisfy 

n 



Lemma 3.2 The multi-index factorial powers of degree \a\, {mh)'^^\ satisfy 

d^^{mh)^^^ ^a\5^^p for\P\ = \a\. 

Lemma 3.3 The multi-index factorial powers [mh)'^^^ of degree \a\ approxi- 
mate the classical multi-index powers x*^"-' of degree \a\ at each point x = rah. 

Let us remark that we have the following relationships between the multi- 
index factorial powers and the usual powers: 

Theorem 3.1 The powers (rnih)°'' and {rriih)'^' are related by 

OLi 

where S^^ are the Stirling numbers of the first kind and T^.' are the Stirling 
numbers of the second kind. 

The sketch of the proof of this theorem can be found, e.g., in [TT] . 
Theorem 3.2 The multi-index powers {mh)^°-^ and {mh)'"^^ are related by 

a 

{mhf^^ = ^ K'^p{mhiP\ (6) 

|/3|=0 

(™/,)(") = L'^3{mh)f. (7) 

|/3|=0 



Moreover, 



l/3| h 

^t- E E ••■E^r^r^.-'-^r^i-z..,- 

l/3| i„-i b 

P / J / J / J 'l '2 — 'l |P|— 

/„_i=o;„_2=o (1=0 
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We will just prove identity (EJ. The proof of identity is analogous to the 
proof of identity . 

Proof: Using Theorem 13.11 and multiplying the polynomials (mi/i)"^ and 
(7712/1)"^, we obtain 

ai+a2 
/3i+/32=0 

with R'/5l^/32 _ CQI CQ2 A 

Using again Theorem 13. II and multiplying the polynomials {mih)°'^ (7712/1)"^ 
and (7713/1)"^, we obtain 

ai+c(2+a3 

(mi/l)"i(m2M"^(™3/i)"^ = ^ ^fj;f^2:a3('™l^)f^("^2MT'^("^3/»)^^'^ 

/3i+/32+/33=0 

Wltil -f^ai, 02,03 ~ ^/2=0 2^h=0^h ^l2-llfJl+p2+f33-l2- 

Applying this procedure recursively, we obtain 



(m/.)(") = ^X^(m/z)!p"' (8) 

|/3|=0 

With if," = e1!L,=o e t :=o • • • e to sr:_,^ ■ ■ ■ s^u^,- ■ 

For all what follows, let 11^ denote the space of all Clifford- valued polyno- 
mials of degree k, P^, generated by the powers (m/i)^'' of degree \a\ — k, and 
IT'= be the countable union of all Clifford- valued polynomials of degree fc > 0. 
Furthermore, let Aif — n kerD^ be the space of discrete monogenic poly- 
nomials of degree fc. Based on Lemma 13.11 13.21 and 13.31 we will show that it 
is possible to obtain discrete versions for the Fischer decomposition as well as 
define discrete versions of the Euler and Gamma operators. 

3.2 The main theorem 

For two Clifford- valued polynomials of degree fc, and e n± given by 

\a\—k 
\a\—k 

we define the Fischer inner product by 

[Pk^QkU--^ E «!Sc(^5±). (9) 

\a\=k 
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Denoting by -P^(-D^) the difference operator obtained from the polynomial 
in powers of mh by replacing rriih by the difference operator dj^^ (just like 
in the continuous case, c.f. j2]), we have by Lemma l3?2l the identity 

[P±, Q±], Sc (lfiDf)Q^m Pt, e n±. (10) 

With other words, we can express the Fischer inner product by applying the 
difference operator P^{D^) to the polynomial P^ and evaluate the scalar part 
at the point rah = 0. 

Moreover, due to — —D^ the Fischer inner product has the important 
property: 

[{mh)Pt, Q^U = -[Pk, D^QkW (11) 
This property allows us to prove the following theorem: 

Theorem 3.3 We have 

= Mf + {mh)Ut,- 

Moreover, the subspaces Ai^ and {mh)Il^_^ are orthogonal with respect to the 
Fischer inner product. 

Before we prove Theorem 81 we will prove the following inclusion property: 
Lemma 3.4 For the set D^^^ ■= {D^Pk^ ■ Pk^ G Hfc I' '^^ have the inclusion: 

^?n± := {D^p^ : p± e n±} c n±_,. 

Proof: 

Let P^{mh) — J2\a\=ki''^^)T ^ "-a £ Ilfc • Applying -D^, we obtain from 
9^*(m/i)^^ — ai{mh)'^ the identity 

n 

{D^P^){rnh) ^Y.T. ("^M^"""^«.e,a± (12) 

Because aiCia^ is a Clifford constant we have a linear combination of poly- 
nomials of degree |a — 8^1 = /c — 1 on the right hand side of Ijl2|l . Hence, 

Proof: (Theorem ESI Because of = {■mh)n^_^ + ((mft.)n^_ it is 
enough to prove that ((m/i)n^_j^)^ — A4f_^. For this we choose P^_i G ^f-i 
arbitrarily and assume that for some P^ £ 11^ we have 

[(m/i)P±_„P±]^-0. 

Due to ® we have [Pk-i,D^P^]h = for all P^_^. As D^P^ e U^^^ by 
LemmaEH we obtain D^P^ = or P^ e M"^ . This means that ((m/i)n^_j)^ C 
Mk- Now, let P^ e M^. Then we have for each P^_^ e n^_^ 

[{mh)Pt„Pt]H = -[Pti^D^P^U 
= 
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and, therefore, ((m/i)n^_^)^ — A^^_i. ■ 

From this theorem we obtain the Fischer decomposition with respect to our 
difference Dirac operators . 

Theorem 3.4 Fischer decomposition Let G IlJ then 

k-l 

P^imh) = J2i^ihyMtA^h). (13) 

s=0 

where each denotes the homogeneous discrete monogenic polynomials of 
degree j with respect to the Dirac operators . 

3.3 DifFerence Euler and Gamma operators 

Based on Lemma l^TTl we will introduce discrete versions of the Euler and Gamma 
operators presented in 0. 

First of all, we introduce the second order difference operator by 

n 

Ai-ThJ2{rn,h)d^'dr. (14) 

i=l 

Definition 3.1 For a lattice function fh '■ K/J; CTo.n, we introduce the differ- 
ence Euler operator by 

n 

{E^fh){mh) = ^(m,/i)(5±7/,)(m/i T he,) 

i=l 

and the difference Gamma operator by 

{T^fh){mh) = -Y,eML%fh){mh) - {A^fu){mh), 
j<k 

where L% {m,h)di^ - {mkh)dfj . 

It looks surprising that we have in the definition of the Gamma operator 
a term which contains second order differences, but we would like to remark 
that for ft, ^ this term vanishes and we will get the usual continuous Gamma 
operator. As a matter of fact this term arises due to the fact that in the 
discrete case translations are involved in the definition of finite differences/ 
finite difference operators. 

Using the definition of the difference Euler operator and Lemma 13.11 we 
obtain for polynomials homogeneous of degree k, P^ e Tlf, E^P^ — kP^ , 
and, moreover, we can show that a function fh homogeneous of degree k satisfy 

fh = kfh- This fact provides a good motivation for calling E'^ Euler opera- 
tor, i.e. an operator who measures the degree of homogeneity of a homogeneous 
function. 
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It follows from the definition of the Euler and Gamma operator that 

n 

{mh)Dih - -Y,{m.h)d^'h + Y.e,ekL%h (15) 

1=1 j<k 

= -iE^+Tt)h (16) 
The proof of ifTfijl is easily obtained by adding and subtracting the term 

n 

-^{mih)d^\fh{mh) - fh{mh ^ hei)), 

i=l 

which is the same as fh, on the right hand side of Ijl5|l . Moreover, for discrete 
monogenic polynomials of degree k, e -^fc i we have F^Af^ = —kM^. 
For all what follows, we introduce the difference operators 

n 

B^ = ±hY,9^\ (17) 

n 

Chfh = Y.^m,h)e,}h{- T he,) (18) 



i=l 

Rtr-^I + E^-^t (19) 

VH%^Rtr + lBt (20) 

where / is the identity operator and r a real number. 
From the identity 

n 

i=i 



we get 



"2-R^,„/2/'i 



{D^imh)) h - (-2i?tv2 + + r^A, (21) 

by applying identity (jl(ill . 

Now, we will show some important facts regarding our difference operators. 



Proposition 3.1 For a lattice function fh ■ K|J C£o,n, we have 
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Proof: 

Starting from the definition, we can split D^Ej^ fh in the sum 

{Df^E^h){mh) = I^{mh)+lt(mh) (22) 

with 

n 

and 

n 

Applying the product rule for finite differences ^ in /f*^ , we obtain 



Itimh) = ^ e, {{d^'h){mh) + {m,h){d^^^\fh){mh T he,)) 

n 
i=l 

On the other hand, 

n 

■-(mh) = ^ e, J2{m,h) [{d^'d^' h)(mh T /le,)) . 



1^2 



Thus, we have 



{D^Eih){mh) = {Dfjh){mh) + ^ e,(m,/i)(9±*a,^^' A)(m/i T /le,) 
= (Z?±/,0(m/i) + (i?±Z?±/,0(m/^). 



Proposition 3.2 For a lattice function fh ■ KJJ C£o,n, we have 

D^{{mh)h) - -2F±„/2A - {mh)Dih (23) 

Proof: Using the product rule for finite differences (j2l and the identity 

—2mih = ei{nih) + (jnh)ei, z = 1, . . . , n, we get 

n n 

i?±((TO/i)A(m/i)) = -^A(m/i±/ieO-2^(m,/i)a±7,.(m/i) 
-(m/i)(i^±/,0(m/i) 

n 

= -nfhimh) - ^{fh{mh ± hei) - ,fh{mh)) - 
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-2(E±A - Aj^h){mh) ~ {mh){D^h){mh) 
= -nh{mh) - 2{E^h " A^h){mh) - 
-{Bfjh){inh) ~ {mh){D^h)imh) 

= -2(i?±„/2/, + ^B^h){mh) - {mh){D^h){mh) 

= -2(Vt„/2/'0("^^) - {rnhKD^fk){mh). 

m 

From Pronosit,ion l3 . 2l a,nd from the commutation properties D^A^ ~ ^h^)^ 
and D^^B^ = B^D^ follow the operator relations 

D^Rtr-Kr+lDt (24) 

Combining Proposition with operator relation H25II . we have for M^_^ S 
Mts, 

= i-^fvL/2Mts- (26) 

and 

{D^nimhfMtj = {Df^f(-2V^^^^,{{mhfMt,) 

= i-^)X%2+2\%2Mts 

+2iDif {{mh)V;^,^^,iimhfMts))) 

~ V ^) y^h,n/2+2^h,n/2 ^ h,n/2+l^ h,n/2 

= i-'^fyL/2+iyL/2^ik-s m 

Continuing this procedure, we obtain by recursion 

{D^Y^imhYMt,) = (-2)^T/±„/2+,_2C/2+.-3---<„/2A^t.- (28) 

From this follows also [mhy M^_^ G ker(_D^)*+^. Formula l|28|l gives us 
a motivation to find expHcit formulae for the polynomials M^_^. To this end 
we need an explicit formula for the inverse of the iterated composite operator 
V.^ /n I tV.^ ,„ , , . . . V,^ This means that we have to find an explicit 

/i,n/2+s — 2 n,n/2+s — 3 h,n/2 ^ 

formula for the inverse of the operator V^^.. Unfortunately, we are only able to 
get an expHcit formula for the operator R^^ (c.f. 0). 
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Theorem 3.5 For a lattice function fh ■ KJJ C£a,n and for r > 0, the 
difference operator J^^ defined by 

iJtrMimh)= J2 hdi({thTh)%-^^fH{{th){mh)) 
t/ie[o,i]h 



satisfies 

Hereby we denote [0, 1]^ — [0, [0, 1]^ = (0, and 



Proof: (c.f. [9]) For fh : MJJ ^ C4,« and r > 0, 

fh{mh)^ hd^({th^h)%^fh{{th){mh)) 
the[o,i]h 

By a direct calculation, 

4 {{thTh)%^fh{{th){mh))) 

= r{th T h)%-^^ fh{{th){mh)) + {th)%\dfjh){{th){mh)) 
= {th T h)%~^^ {rfh{{th){mh)) + th{d^fh){{th){mh))) 

On the other hand, applying the difference version of the chain rule and the 
relation Y^^=ii^''^i^)'^h^ — ~ ' '^^ obtain 



th{d^fh){m{mh)) = Y.{th){m,h){d^^ fh){{th){mh)) 

1=1 

= {E^Jh)m){mh)) - {AfMmimh)) 

Therefore, 

fhimh) = r{J^Jh)irnh) + {E^J^Jh){mh)~{Aijijh){mh) 
= (Rtr^trMimh). 

From the above two identities and the definitions of R^^ and J^^, we get 

m 

Now, the construction of the inverse for Vj^^ seems to be obvious. But, the 
obvious choice 

the[o,i]± 
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is not an inverse of Vf^^., which can be easily checked in the following way. 
If we use the same technique as above, we obtain 

h{mh)= hd^ {{th)^r^^ M{th){mh))) 

t/ie[o,i]J 

and by direct calculation 

4 {{th)^:^^hi{th){mh))) 

- r{th)^^-''> fhithx) + [th ± h)%\d^fh){{th){mh)) 

= {th)%-^^ {rh{{th){mh)) + {th ± h){d^fh){{th){mh))) 

On the other hand, applying the difference version of the chain rule and the 
relation J2^=ii''^i^)^h^ ~ " '^^ obtain 

n 

{th±h){d^h){{th){mh)) - ^((i/i)K/i)±/iTO,)(af^7/.)((i/i)(m/i)) 

2 = 1 

= {EfjHmh){mh)) - {Afj^){{th){mh)) 

n 

±hJ2m.Hdi^h)iith){mh)), 

i=l 

but 

n n 

1=1 i=l 

= {BtM(th){mh)). 

3.4 Difference operator calculus 

Now we will establish some properties for our difference operators introduced 
in Section rOI 

Using the difference properties 

and 

{m^h)d^'d^\mh T /le,)^"^ = {a, - l)im,h)d^\mh T he.f^'^ 

we obtain by direct calculation the following formulae for homogeneous polyno- 
mials of degree fc, G 11^, 



S±P± = ±khP^, (29) 
l±h' 



^kPk = T^Pk^ (30) 



kh^ 
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Now, using the difference rules ^ and we get 





- {7nh)B^fh + hl^Umh) + h^D^fh, 


(33) 




= {mh)Cth-E^h, 


(34) 






(35) 


Ei{{mh)fu) 




(36) 




= {mh)Rljh + (1 ± /i)C±A, 


(37) 




= imh)V^rh + (1 ± h)C^h + 
+ ^{hl^fhimh) + h^D^fh), 


(38) 



where 1^ := ±I]r=i 

As a direct consequence of formulae Ij32ll and Ij36ll , we obtain for the discrete 
homogeneous monogenic polynomials of degree k, E M^, 

Tt{{mh)M^) = -E^{{mh)M^)~{mh)D^{{mh)M^) 

by applying Theorem 13.21 and relation ifTfill . 

Moreover, using identities (ED, lEHIl , JSU , ISll , CEl and Proposition 13.21 we 
obtain the relation 

{D^imh)) {{mh)M^) = (-2<„/2 + + r±)((m/.)M±) 

= ±hk{mh)M^ - {2±2h)C^M^. (40) 

Applying relations ll36)l and l(3ijl we have 

E^iimhfM^) = k{mhfMf + 2{mh)C^M^ - kM^ . (41) 

From Theorem 13.21 and formulae l(33jl and , we get 

Df^{{mKfMi) - -(2 ± 2K)CiMi + hl^M^. (42) 

Using Ijl6|l and formulae Ij41|l and Ij42ll . we obtain 

Tf^{{mhfM^) - ~k{mh)Hl^ ± 2(m/i)C±A/± + fcM± - hl^{mh)M^. (43) 

Applying recursively our formulae Il33ll - lj38|l . it is possible to obtain explicit 
formulae for D^{{mhYMf), E^{{mhYMf) and T^{{mhyM^) by induction. 
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3.5 Homogeneous powers 

Contrary to the continuous case, the classical product between the variable mh 
and the homogeneous polynomial is not homogeneous. However, apply- 
ing formulae lO and proved in Theorem 13.21 we can say that the product 
{mh)P^ can be expressed as a linear combination of homogeneous polynomials 
up to degree k -\- \. On the other hand, the powers x'^ = {mhy are not homo- 
geneous. For this purpose, we will introduce the discrete analogues of a;* in the 
following way: 

Starting from the multi-index factorial powers of degree {mhf-^\ we 
introduce the polynomials Hf , s G N by 

i|=fc 



Hf{mh) 



0- 



As a direct consequence of the identity 



(44) 



{-lf\x\'^^ if s ^ 2k, k eNq 
xx^'' i/ s = 2fc + G No 

we can conclude by Lemma and by the multinomial theorem, that the poly- 
nomials Hf^ give rise to homogeneous polynomials of degree s which approxi- 
mate the powers = {mhy for small mesh width h > 0. 

As a direct consequence, the operator formulae proved in Subsection l3.4l are 
fulfilled for the powers and, moreover, by direct calculation, we obtain the 
additional properties 

and 

Let us remark that the term C^Hf is the discrete version of the multiplication 
xx'^ in the continuous case. 



4 A discrete harmonic Fischer decomposition 

According to the classical theory of the finite differences, the usual approxima- 
tion of the Laplacian is given by 

u(mh -\- hei) + u(mh — hei) — 2u(mh) 
{Ahu){mh) = 2^ — 

i=l 
n 

= Y.(9rdt^Kmh). (45) 

i=l 

The first problem that arises now is that not all of our partial difference 
operators do commute in the certain sense (c.f. 'W.) and, moreover, we have 
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no factorization of the discrete Laplacian Ah by means of our difference Dirac 
operators considered above, that is D^Df^ ^ — eoA/j. 

Let us restrict ourselves in this section to the case of quaternion-valued 
functions defined on lattices in . 

Let us remark that the quaternionic variable mh is identified with the 4x4 
matrix 

—rriih 


m^h 
—m2h 



( 



mh 





mi ft, 

7712 ft 

V "73ft 

In (3 for a lattice function fh 

3 



—7772/7 



7771/7 

II given by 



-7773/7 ^ 
7772 ft 

777lft 




VecA 



a finite difference approximation of our Dirac operator was defined in the form 



( 







V 



-dl 



-dl 



-div^ Vec//i 

' curl+VecA 



/ 



fl 
fl 
\ffj 



(46) 



-dl 



dl 

dl 
V dl 

-div^Vec/ft 
gi'ad^/h +curl;^Vcc/,i 



di: 



-dl 

d-^ 





f ^° \ 




fl 




fl 


) 


\fl 



(47) 



with div±VecA = Ell di^'fl, grad±/o 



ELi(5^7.°)e. and 



curl^ Vec fh = 



D.f = 



ei 62 63 
^? ^? 

//i fh fh 

In the latter form one can easily see the similarity with the usual Dirac 
operator 

— divVec / 
gradSc / + curlVec / 

Using the discrete identities 



div^curl^Vec/,i 







curl±grad±/° = 
curljJ=curlJVec/;i = -AhVecfh + gradJdiv^Vee //i 
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we obtain the following factorization of the discrete Laplacian 

7/,. = 



^A,,VecA 



(48) 



Now, we are able to obtain a Fischer decomposition for the discrete Dirac 

+- 

h ■ 



operators Dj^ and D 



Using the fact that 

D7 





jrad^ [mh) 




(a) 



ai{mh) 
a2{mh) 
\ a3{mh) 



(a— ei) 

+ 

(0-62) 

+ 



as well as 



D+-{mh) 



(a) 



jiad^ (m/i) 



(a) 



ai[mn)_ 
a2{mh)'f 



we can prove as in Lemma \l^A\ the inclusion properties D'^ 11^1" C 
Df^^H^ C and, moreover, replacing Df^~ by and D^'^ by I?^ in the 

inner product ifTTijl . we obtain the Fischer decompositions: 

Theorem 4.1 Fischer decomposition for and 

Let PjT G ( respectively, G 11^ ) then 

fe-i 



(49) 



s=0 

E 

s=0 



[mhYMl 



(50) 



where each ( respectively, ) denotes a homogeneous discrete mono- 

genic polynomial of degree i , that is, MJ^ G IVJ r\kcT:: D'j^^ (respectively, G 

n+nkeri?+-;. 

From the factorization property H48II . we have 

[(m/j)2F±,Q±]„ = -[P±,A„Q±],„ 
which allows us to obtain the Fischer decomposition for the discrete Laplacian: 
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Theorem 4.2 Fischer decomposition for 



Let e the 



2s<k 



where each Hj denotes a homogeneous discrete harmonic polynomial of degree 
j, that is, Hf eUf n ker A^. 

As a consequence of Theorem 14. II we obtain Fischer decompositions which 
relate the discrete harmonic and the discrete monogenic polynomials. 

Corollary 4.1 Fischer decomposition Let G 11^ flkerA/i then 



M+- + {mh)M+Z^. 



(51) 
(52) 



where each ( respectively, ) denotes a homogeneous discrete mono- 

genic polynomial of degree j , that is, £ IVJ r\kei: D'f^^ (respectively, AIj € 

n+ n ker 



To define the Euler and Gamma operator , and Ej^ 



,F^+ for the 



modified Dirac operators D^" and D^^ , respectively, we start to calculate 
the products {mh)D'^^ fh and {nih)D^^ fh- By straightforward calculations we 
obtain 



{mh)D-+fh 



a-2 j^O 
"h Jh 



ei 

mih 



62 




eo e2 es 

mih m^h 



eo 

— TO2^ 



ei 

-dl 



63 

mih 



eo 6i eg 
7713/1 mih 



dZ 



fl 



fl 



(53) 



and 



{mh)D+-fh 



( dlfl 
dlfl 
dl'fl 
V dl'fl 



a2 i^O 
"hJh 

'h Jh 

dlfl 
11 



dlfl 
dZ'fl 
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ei 


62 


63 




60 


62 


63 


mih 


■m2h 






mih 




777-2/1 


dl 




dl 




dl 


d-n' 


d-H' 



fh 



+ 



+ 



60 61 63 

-7772/1 —m^h mih 



-dl 



60 61 62 

7773/7 77^2/1 '771/7 



5: 



d: 



fh- 



(54) 



Hence, we can define the difference Euler operators E^^ ^ and E^^ as 

{E^+h){mh) = 
( (dj^'mimh + hei) 

{dlfl){mh-hei) 
{dlfl){mh-her) 



(9,-V^°)(m/7 + /7e2) (a,-V^°)(m/7 + /7e3) \ 



and 



{E+-h){mh) = 

( {dlOimh-hBi) 
{dlfk){mh-hei) 

\ {d^'fl){mh + her) 



{dlfl)imh-he2) 
{d^^mmh + he^) 
{dlfDimh-he^) 



{dlfl){mh-he:,) 
{d-^f^){mh + hes) ) 



idlfjl){mh - he^) 
{d^\fl){mh + he2) 

{dlfl){mh - he,) 
{d^^fl){mh + he2) 



The difTerence Gamma operators and are defined by 

{Tl+h){mh) = 

( {d^'dlfl){mh) {d^^dlfl){mh) 
-{dfdlfl){mh) 



h 



\ 



id^'dlf^){mh) 
idh'dUDimh) 
-(.d^'dlfl){mh) 
-{d^'dlfl){mh) 



-{d^^dlfl){mh) 



id^'dlf,)imh) 
-{d^'dlfl){mh) 
-{d^'dlfl){mh) 
{d^'dlfl){mh) 



61 


62 63 


mih 


7772/7 7773/7 


d-u' 


dn' d-,' 


60 


61 


—7772/1 


—m^h n 


-df 


-dl 


60 


61 62 


m^h 


7772/7 mih 


d-h' 


dl di 



60 62 63 
mih m^h 7772/7 



d-u' 



dl 



dl 



flimh) 



63 



dl 



flimh) 



fl{mh) 



and 



mih 
m^h 
m^h 



{dlO{mh-hB,) 
id-^fl)(mh + he,) 
{d^^fl)(mh + he,) 

{dlfl){mh-he^) 



mih 
7772/7 
m^h 



7774/7 
7772/7. 
m^h 



{r+-h){mh) 
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/ (.d-'dlf,)imh) 



-{d^'dlfl){mh) 



-{d^'dlfl){mh) 
{dT,'dlfl){mh) J 



mih 



= -h 



V -{du'dlfl){mh) 



ei 62 63 
— mih ni2h m^h f^{mh) — 



al a2 a3 




Go ei 63 
— —7772/1 —m^h mih f'^{mh) 



Go ei 62 

— m^h 7772/7 mih ff^{mh). 



"h "h 



As in Subsection we have {mh)Dj^^ = —E^^^ ^'^h^ (respectively, 
{mh)D+- = - r+") and the polynomials e satisfy E^+P^ = 

kPjT , (respectively, E^^ P^ — kPj^). Moreover, if PjT G kerD,""'' (respec- 
tively, Pft G ker_D^^) then we have T^+P^T = —kPf^, (respectively, rjJ^^Pj^ = 



Like in Theorem 13.21 we can prove the operator property D'f^^Ef^^ ~ I + 
E^^ Dy^^ (respectively, D'^^ E^^ = I + E'^^ D^^). In the same way we get 
analogous relations to the ones presented in Subsection l3.3l and in Subsection l3.4l 
In addition it is also possible define the discrete versions of the quaternionic 
powers {mhy with respect to our difference Dirac operators DJ^^ and D'^~ , 
using a similar construction as in Subsection 13.51 

At this point it would be interesting to know if the operator setting we 
discussed here for the quaternionic case has an equivalent operator setting in 
the general case of Clifford algebras. Up to know it is not known, but we will 
discuss it in a forthcoming paper 0]. 
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